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. Abstract. We present the first gas dynamical simulations of the photoevaporation of cosmological 

pLn ' minihalos overtaken by the ionization fronts which swept through the IGM during reionization in 

a ACDM universe, including the effects of radiative transfer. We demonstrate the phenomenon of 
I-front trapping inside minihalos, in which the weak, R-type fronts which traveled supersonically 
across the IGM decelerated when they encountered the dense, neutral gas inside minihalos, becom- 
ing D-type I-fronts, preceded by shock waves. For a minihalo with virial temperature Ty[^ < lO^K, 
^ ■ the I-front gradually burned its way through the minihalo which trapped it, removing all of its bary- 

Q> I onic gas by causing a supersonic, evaporative wind to blow backwards into the IGM, away from the 

i exposed layers of minihalo gas just behind the advancing I-front. 

CO ' Such hitherto neglected feedback effects were widespread during reionization. N-body Simula- 

tions and analytical estimates of halo formation suggest that sub-kpc minihalos such as these, with 
^—^ . Tvir < lO^K, were so common as to dominate the absorption of ionizing photons. This means that 

previous estimates of the number of ionizing photons per H atom required to complete reionization 
which neglected this effect may be too low. Regardless of their effect on the progress of reionization, 
• however, the minihalos were so abundant that random lines of sight thru the high-z universe should 

Qh' encounter many of them, which suggests that it may be possible to observe the processes described 

Q . here in the absorption spectra of distant sources. 
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>< : I-FRONTS AND MINIHALOS AT HIGH REDSHIFT 

The first sources of ionizing radiation to condense out of the dark, neutral, opaque IGM 
reheated and reionized it between z ~ 30 and z ~ 6. Weak, R-type ionization fronts sur- 
rounding each source swept outward thru the IGM, overtaking smaller-mass virialized 
halos, called minihalos, and photoevaporating them. High-redshift sources of ionizing 
photons may have found their sky covered by these minihalos. The photoevaporation 
of minihalos may therefore have dominated consumption of ionizing photons during 
reionization! 1 ]. In this paper, we focus on the currently favored ACDM model (Qq = 0.3, 
Ao = 0.7, h = 0.7, Q.hh^ = 0.02). In this model, the universe at z ^ 6 was already filled 
with minihalos capable of trapping a piece of the global, intergalactic I-fronts, photoe- 
vaporating their gaseous baryons back into the IGM. Prior to their encounter with these 
I-fronts, minihalos with Tvir < lO^K were neutral and optically thick to hydrogen ion- 
izing radiation. N-body simulations in a cubic volume lOOkpc on a side (see Fig. la,b) 
reveal that at z = 9 minihalos with T^^^ < lO'^K (i.e. M < lO^-^M©) were separated on av- 



erage by only d ~ 7kpc while their geometric cross section together covered ~ 16% of 
the area along every 100 kpc of an average line of sight [5]. If the sources of reionization 
were large mass halos with T > Tyir, then these were well-enough separated that typical 
reionization photons were likely to have been absorbed by intervening minihalos. 

To demonstrate this in a statistically meaningful way with more dynamic range than 
N-body results, we combine the Press-Schechter (PS) prescription for deriving the aver- 
age number density of halos with the truncated isothermal sphere (TIS) halo model[2,Q], 
to determine which halos are subject to photoevaporation and how common they are, as 
function of their redshifts. We then compute the fraction of the sky i^cover,source5 as seen 
by a source halo of a given mass, which is covered by opaque minihalos located within 
the mean volume per source halo. If halos with M ^ lO^M© are the reionization sources, 
their minihalo covering fraction is close to unity and increases by a factor of a few if we 
take account of the statistical bias by which minihalos tend to cluster around the source 
halos (see Fig. lc,d). 



100 kpc box, TviR > 10*K 100 kpc box, Tv,r < 10*K 




FIGURE 1. Dark matter halos of mass M > IO'^-^Mq (i.e. Tyi^ > lO'^K) (a) (top left) and IO^ ^^Mq > 
M > lO^-^M© (i.e. Tvir < lO'^K) (b) (top center) found in P^M simulation of ACDM universe in 100 kpc 
(proper) box, at z = 9. (c) (top right) Fraction of sky covered by minihalos located within the mean 
volume per source halo versus halo mass, at different redshifts, computed using the TIS model and 
the PS approximation, corrected for linear bias, (d) (bottom) Sky covering by minihalos as seen from 
a 1.1 X IO^Mq source, from P^M simulation with 50 kpc box at z = 9. All minihalos within 25kpc are 
plotted. The covering fraction is 23.6%. 



THE PHOTOEVAPORATION OF MINIHALOS 



We have performed radiation-hydrodynamical simulations of the photoevaporation of a 
cosmological minihalo overrun by a vv'eak, R-type I-front in the surrounding IGM, cre- 
ated by an external source of radiation[3]. Minihalos are modeled as TIS of CDM + 
baryons We consider 3 different source spectra: (1) QSO-like: Fy °^ v ^'^ (v > 

Vh); (2) Stellar blackbody: Teff = 50, 000 K; (3) "No Metals" Stellar T^ff = 100, 000 K. 
The I-front encounters the minihalo at redshifts zinitiai = 

(6,7,9,11). The flux lev- 
els are Fq = %,56(v > V//)/r2jp^ = (0.1,0.5, 1,2,5, 10, 1000) (where % = A^ph,56 x 

10^^ ionizing photons cm^^s"', and r = rMpc x IMpc). The halo masses are M^aio = 
(10^, 10^ 10^ 10^2 X 10^4 X 1O^)M0 (where 4 x lO^M© corresponds to T^i, = \0^K 
at z — 9). The minihalo shields itself against ionizing photons, traps the R-type I-front 
which enters the halo, causing it to decelerate inside the halo to close to the sound speed 
of the ionized gas and transforms itself into a D-type front, preceded by a shock. The side 
facing the source expels a supersonic wind backward toward the source, which shocks 
the IGM outside the minihalo. The wind grows more isotropic with time as the remain- 
ing neutral halo material is photoevaporated. Since the gas itself was initially bound to a 
dark halo with Oy < lOkms^^, photoevaporation proceeds unimpeded by gravity. Fig- 
ure 2 shows the structure of the photoevaporation flow 60 Myr after the global I-front 
first overtakes a 10^ Mq minihalo. 

The evaporation time per halo, t^^, and the number, , of ionizing photons absorbed 



t = 60 Myr, BB 5e4, (lO^M^.Q,!), (0.3.0.7,0.7), (25,3), w/l-front 




FIGURE 2. One time-slice, 60 Myr after I-front caused by source (located far to the left, along the x- 
axis) overtakes a lO^M© minihalo [centered at (r,x) = (0, 1 .06 x lO'^'^cm)] at z = 9 in the ACDM universe, 
for source with Fq = 1, a stellar BB spectrum ~ 50, 000 K. (a) (top) Isocontours of atomic density, 
logarithmically spaced; (b) (bottom) Velocity arrows are plotted with length proportional to gas velocity. 
An arrow of length equal to the spacing between arrows has velocity 25kms^^; minimum velocities 
plotted are 3kms^'. Solid line shows current extent of gas initially inside minihalos at z = 9. Dashed line 
is I-front (50% H-ionization contour). 



per minihalo H atom during this evaporation time have important implications for 
the reionization of the universe. We summarize in Figure 3 the dependence of these 
quantities on the value of the Mhaio, Zinitiah and the source spectrum, based on our 
simulation results. 
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FIGURE 3. Photoevaporation times ?ev for individual minihalos (top panels) and total number of 
ionizing photons absorbed per minihalo H atom during this photoevaporation (bottom panels). Different 
panels show the variation in these quantities with the input parameters which label the horizontal axes: 
Miaio = minihalo mass, z = redshift when I-front first encounters minihalo, and Fq = dimensionless 
ionizing photon flux, with spectrum which is either a stellar blackbody, with = 50,000K, (circles) 
or Jh. = 100,000K (squares), or else QSO-Uke (triangles). 



